Mathematical Modeling of ceRNA-Based Interactions.

Competing endogenous RNA (ceRNA) molecules have emerged as key players in regulating gene expression, increasing the complexity of the range of possible dynamics within a cell. The actions of competing RNA typically are sponging behaviors, in a manner that fine-tunes gene expression, but there are particular network structures that may show destabilization due to ceRNA interactions. In this chapter, we discuss how these interactions can be modeled and probed from a mathematical, first-principles perspective.

Pseudogenes as Competitive Endogenous RNAs: Target Prediction and Validation.

Pseudogenes may regulate expression of their parental genes as well as other protein-coding genes through various mechanisms. One such mechanism is the ability to act as competitive endogenous RNA (ceRNA) and participate in microRNA-mediated cross-regulation. Here, we outline how to predict the targets of pseudogene ceRNAs bioinformatically and how to validate them experimentally.

Pseudogenes as Competitive Endogenous RNAs: Testing miRNA Dependency.

Pseudogenes may function as competitive endogenous RNAs (ceRNAs), where they regulate the expression of genes by sequestering shared miRNAs. ceRNAs are becoming more extensively identified and studied, and demonstrating the dependence of their effects on miRNA sequestration is critical to establish them as ceRNAs. Here, we outline an experimental approach to assess the miRNA dependency of a candidate pseudogene ceRNA.

An Overview of the Computational Models Dealing with the Regulatory ceRNA Mechanism and ceRNA Deregulation in Cancer.

Pools of RNA molecules can act as competing endogenous RNAs (ceRNAs) and indirectly alter their expression levels by competitively binding shared microRNAs. This ceRNA cross talk yields an additional posttranscriptional regulatory layer, which plays key roles in both physiological and pathological processes. MicroRNAs can act as decoys by binding multiple RNAs, as well as RNAs can act as ceRNAs by competing for binding multiple microRNAs, leading to many cross talk interactions that could favor significant large-scale effects in spite of the weakness of single interactions. Identifying and studying these extended ceRNA interaction networks could provide a global view of the fine-tuning gene regulation in a wide range of biological processes and tumor progressions. In this chapter, we review current progress of predicting ceRNA cross talk, by summarizing the most up-to-date databases, which collect computationally predicted and/or experimentally validated miRNA-target and ceRNA-ceRNA interactions, as well as the widespread computational methods for discovering and modeling possible evidences of ceRNA-ceRNA interaction networks. These methods can be grouped in two categories: statistics-based methods exploit multivariate analysis to build ceRNA networks, by considering the miRNA expression levels when evaluating miRNA sponging relationships; mathematical methods build deterministic or stochastic models to analyze and predict the behavior of ceRNA cross talk.

LuxT controls specific quorum-sensing-regulated behaviors in Vibrionaceae spp. via repression of qrr1, encoding a small regulatory RNA.

Quorum sensing (QS) is a process of chemical communication bacteria use to transition between individual and collective behaviors. QS depends on the production, release, and synchronous response to signaling molecules called autoinducers (AIs). The marine bacterium Vibrio harveyi monitors AIs using a signal transduction pathway that relies on five small regulatory RNAs (called Qrr1-5) that post-transcriptionally control target genes. Curiously, the small RNAs largely function redundantly making it difficult to understand the necessity for five of them. Here, we identify LuxT as a transcriptional repressor of qrr1. LuxT does not regulate qrr2-5, demonstrating that qrr genes can be independently controlled to drive unique downstream QS gene expression patterns. LuxT reinforces its control over the same genes it regulates indirectly via repression of qrr1, through a second transcriptional control mechanism. Genes dually regulated by LuxT specify public goods including an aerolysin-type pore-forming toxin. Phylogenetic analyses reveal that LuxT is conserved among Vibrionaceae and sequence comparisons predict that LuxT represses qrr1 in additional species. The present findings reveal that the QS regulatory RNAs can carry out both shared and unique functions to endow bacteria with plasticity in their output behaviors.

Massively Parallel Analysis of Regulatory RNA Sequences.

The stability of RNA transcripts is regulated by signals within their sequences, but the identity of those signals still remain elusive in many biological systems. Recently introduced massively parallel tools for the analysis of regulatory RNA sequences provide the ability to detect functional cis-regulatory sequences of post-transcriptional RNA regulation at a much larger scale and resolution than before. Their application formulates the underlying sequence-based rules and predicts the impact of genetic variations. Here, we describe the application of UTR-Seq, as a strategy to uncover cis-regulatory signals of RNA stability during early zebrafish embryogenesis. The method combines massively parallel reporter assays (MPRA) with computational regression models. It surveys the effect of tens of thousands of regulatory sequences on RNA stability and analyzes the results via regression models to identify sequence signals that impact RNA stability and to predict the in vivo effect of sequence variations.

ACE2 polymorphisms as potential players in COVID-19 outcome.

The clinical condition COVID-19, caused by SARS-CoV-2, was declared a pandemic by the WHO in March 2020. Currently, there are more than 5 million cases worldwide, and the pandemic has increased exponentially in many countries, with different incidences and death rates among regions/ethnicities and, intriguingly, between sexes. In addition to the many factors that can influence these discrepancies, we suggest a biological aspect, the genetic variation at the viral S protein receptor in human cells, ACE2 (angiotensin I-converting enzyme 2), which may contribute to the worse clinical outcome in males and in some regions worldwide. We performed exomics analysis in native and admixed South American populations, and we also conducted in silico genomics databank investigations in populations from other continents. Interestingly, at least ten polymorphisms in coding, noncoding and regulatory sites were found that can shed light on this issue and offer a plausible biological explanation for these epidemiological differences. In conclusion, there are ACE2 polymorphisms that could influence epidemiological discrepancies observed among ancestry and, moreover, between sexes.

Effects of Coronavirus Persistence on the Genome Structure and Subsequent Gene Expression, Pathogenicity and Adaptation Capability.

Coronaviruses are able to establish persistence. However, how coronaviruses react to persistence and whether the selected viruses have altered their characteristics remain unclear. In this study, we found that the persistent infection of bovine coronavirus (BCoV), which is in the same genus as SARS-COV-2, led to alterations of genome structure, attenuation of gene expression, and the synthesis of subgenomic mRNA (sgmRNA) with a previously unidentified pattern. Subsequent analyses revealed that the altered genome structures were associated with the attenuation of gene expression. In addition, the genome structure at the 5' terminus and the cellular environment during the persistence were responsible for the sgmRNA synthesis, solving the previously unanswered question regarding the selection of transcription regulatory sequence for synthesis of BCoV sgmRNA 12.7. Although the BCoV variants (BCoV-p95) selected under the persistence replicated efficiently in cells without persistent infection, its pathogenicity was still lower than that of wild-type (wt) BCoV. Furthermore, in comparison with wt BCoV, the variant BCoV-p95 was not able to efficiently adapt to the challenges of alternative environments, suggesting wt BCoV is genetically robust. We anticipate that the findings derived from this fundamental research can contribute to the disease control and treatments against coronavirus infection including SARS-CoV-2.

Early termination of the Shiga toxin transcript generates a regulatory small RNA.

Enterohemorrhagic Escherichia coli is a significant human pathogen that causes disease ranging from hemorrhagic colitis to hemolytic uremic syndrome. The latter can lead to potentially fatal renal failure and is caused by the release of Shiga toxins that are encoded within lambdoid bacteriophages. The toxins are encoded within the late transcript of the phage and are regulated by antitermination of the PR' late promoter during lytic induction of the phage. During lysogeny, the late transcript is prematurely terminated at tR' immediately downstream of PR', generating a short RNA that is a byproduct of antitermination regulation. We demonstrate that this short transcript binds the small RNA chaperone Hfq, and is processed into a stable 74-nt regulatory small RNA that we have termed StxS. StxS represses expression of Shiga toxin 1 under lysogenic conditions through direct interactions with the stx1AB transcript. StxS acts in trans to activate expression of the general stress response sigma factor, RpoS, through direct interactions with an activating seed sequence within the 5' UTR. Activation of RpoS promotes high cell density growth under nutrient-limiting conditions. Many phages utilize antitermination to regulate the lytic/lysogenic switch and our results demonstrate that short RNAs generated as a byproduct of this regulation can acquire regulatory small RNA functions that modulate host fitness.

Characterization of a female germline and early zygote promoter from the transcription factor bZip1 in the dengue mosquito Aedes aegypti.

BACKGROUND: The wide distribution of Aedes aegypti, the main vector of dengue and yellow fever viruses, currently puts three billion people in the world at risk of infection with these viruses. Continuous transmission of these and other viruses despite aggressive efforts to prevent this emphasizes the need to develop new control strategies. Proposals to control disease transmission based on vector engineering, including both population suppression and population replacement, rely on the development of transgenes under the control of regulatory elements able to drive molecules in a specific tissue, time and strength. METHODS: Here we report the characterization of a promoter active in both the female germline and early zygote, derived from the transcription factor bZip1 in the mosquito Ae. aegypti, using transposon-based methods and RT-qPCR. RESULTS: We generated seven transgenic lines carrying AabZip1-reporter constructs and observed expression in both the ovary and early embryo. RT-qPCR analysis was performed to evaluate transcript expression patterns for each line, confirming that transgenic expression from the AabZip1 promoter largely recapitulated the endogenous expression pattern, albeit the strength of maternal expression appeared to be strongly influenced by chromosomal position. CONCLUSIONS: This study provides a new regulatory sequence that can be useful for generating transgenic lines that can become a tool in vector control strategies.

Reading between the Lines: Utilizing RNA-Seq Data for Global Analysis of sRNAs in Staphylococcus aureus.

Regulatory small RNAs (sRNAs) are known to play important roles in the Gram-positive bacterial pathogen Staphylococcus aureus; however, their existence is often overlooked, primarily because sRNA genes are absent from genome annotation files. Consequently, transcriptome sequencing (RNA-Seq)-based experimental approaches, performed using standard genome annotation files as a reference, have likely overlooked data for sRNAs. Previously, we created an updated S. aureus genome annotation file, which included annotations for 303 known sRNAs in USA300. Here, we utilized this updated reference file to reexamine publicly available RNA-Seq data sets in an attempt to recover lost information on sRNA expression, stability, and potential to encode peptides. First, we used transcriptomic data from 22 studies to identify how the expression of 303 sRNAs changed under 64 different experimental conditions. Next, we used RNA-Seq data from an RNA stability assay to identify highly stable/unstable sRNAs. We went on to reanalyze a ribosome profiling (Ribo-seq) data set to identify sRNAs that have the potential to encode peptides and to experimentally confirm the presence of three of these peptides in the USA300 background. Interestingly, one of these sRNAs/peptides, encoded at the tsr37 locus, influences the ability of S. aureus cells to autoaggregate. Finally, we reexamined two recently published in vivo RNA-Seq data sets, from the cystic fibrosis (CF) lung and a murine vaginal colonization study, and identified 29 sRNAs that may play a role in vivo Collectively, these results can help inform future studies of these important regulatory elements in S. aureus and highlight the need for ongoing curating and updating of genome annotation files.IMPORTANCE Regulatory small RNAs (sRNAs) are a class of RNA molecules that are produced in bacterial cells but that typically do not encode proteins. Instead, they perform a variety of critical functions within the cell as RNA. Most bacterial genomes do not include annotations for sRNA genes, and any type of analysis that is performed using a bacterial genome as a reference will therefore overlook data for sRNAs. In this study, we reexamined hundreds of previously generated S. aureus RNA-Seq data sets and reanalyzed them to generate data for sRNAs. To do so, we utilized an updated S. aureus genome annotation file, previously generated by our group, which contains annotations for 303 sRNAs. The data generated (which were previously discarded) shed new light on sRNAs in S. aureus, most of which are unstudied, and highlight certain sRNAs that are likely to play important roles in the cell.

Interplay between genetics and epigenetics in osteoarthritis.

Research into the molecular genetics of osteoarthritis (OA) has been substantially bolstered in the past few years by the implementation of powerful genome-wide scans that have revealed a large number of novel risk loci associated with the disease. This refreshing wave of discovery has occurred concurrently with epigenetic studies of joint tissues that have examined DNA methylation, histone modifications and regulatory RNAs. These epigenetic analyses have involved investigations of joint development, homeostasis and disease and have used both human samples and animal models. What has become apparent from a comparison of these two complementary approaches is that many OA genetic risk signals interact with, map to or correlate with epigenetic mediators. This discovery implies that epigenetic mechanisms, and their effect on gene expression, are a major conduit through which OA genetic risk polymorphisms exert their functional effects. This observation is particularly exciting as it provides mechanistic insight into OA susceptibility. Furthermore, this knowledge reveals avenues for attenuating the negative effect of risk-conferring alleles by exposing the epigenome as an exploitable target for therapeutic intervention in OA.

HOX gene cluster (de)regulation in brain: from neurodevelopment to malignant glial tumours.

HOX genes encode a family of evolutionarily conserved homeodomain transcription factors that are crucial both during development and adult life. In humans, 39 HOX genes are arranged in four clusters (HOXA, B, C, and D) in chromosomes 7, 17, 12, and 2, respectively. During embryonic development, particular epigenetic states accompany their expression along the anterior-posterior body axis. This tightly regulated temporal-spatial expression pattern reflects their relative chromosomal localization, and is critical for normal embryonic brain development when HOX genes are mainly expressed in the hindbrain and mostly absent in the forebrain region. Epigenetic marks, mostly polycomb-associated, are dynamically regulated at HOX loci and regulatory regions to ensure the finely tuned HOX activation and repression, highlighting a crucial epigenetic plasticity necessary for homeostatic development. HOX genes are essentially absent in healthy adult brain, whereas they are detected in malignant brain tumours, namely gliomas, where HOX genes display critical roles by regulating several hallmarks of cancer. Here, we review the major mechanisms involved in HOX genes (de)regulation in the brain, from embryonic to adult stages, in physiological and oncologic conditions. We focus particularly on the emerging causes of HOX gene deregulation in glioma, as well as on their functional and clinical implications.

Non-coding RNA Associated Competitive Endogenous RNA Regulatory Network: Novel Therapeutic Approach in Liver Fibrosis.

Liver fibrosis or scarring is the most common pathological feature caused by chronic liver injury, and is widely considered one of the primary causes of morbidity and mortality. It is primarily characterised by hepatic stellate cells (HSC) activation and excessive extracellular matrix (ECM) protein deposition. Overwhelming evidence suggests that the dysregulation of several noncoding RNAs (ncRNAs), mainly long non-coding RNAs (lncRNAs), microRNAs (miRNAs) and circular RNAs (circRNAs) contributes to the activation of HSC and progression of liver fibrosis. These ncRNAs not only bind to their target genes for the development and regression of liver fibrosis but also act as competing endogenous RNAs (ceRNAs) by sponging with miRNAs to form signaling cascades. Among these signaling cascades, lncRNA-miRNA-mRNA and circRNA-miRNA-mRNA are critical modulators for the initiation, progression, and regression of liver fibrosis. Thus, targeting these interacting ncRNA cascades can serve as a novel and potential therapeutic target for inhibition of HSC activation and prevention and regression of liver fibrosis.

The changing paradigm of intron retention: regulation, ramifications and recipes.

Intron retention (IR) is a form of alternative splicing that has long been neglected in mammalian systems although it has been studied for decades in non-mammalian species such as plants, fungi, insects and viruses. It was generally assumed that mis-splicing, leading to the retention of introns, would have no physiological consequence other than reducing gene expression by nonsense-mediated decay. Relatively recent landmark discoveries have highlighted the pivotal role that IR serves in normal and disease-related human biology. Significant technical hurdles have been overcome, thereby enabling the robust detection and quantification of IR. Still, relatively little is known about the cis- and trans-acting modulators controlling this phenomenon. The fate of an intron to be, or not to be, retained in the mature transcript is the direct result of the influence exerted by numerous intrinsic and extrinsic factors at multiple levels of regulation. These factors have altered current biological paradigms and provided unexpected insights into the transcriptional landscape. In this review, we discuss the regulators of IR and methods to identify them. Our focus is primarily on mammals, however, we broaden the scope to non-mammalian organisms in which IR has been shown to be biologically relevant.

Identification of a Xist silencing domain by Tiling CRISPR.

Despite essential roles played by long noncoding RNAs (lncRNAs) in development and disease, methods to determine lncRNA cis-elements are lacking. Here, we developed a screening method named "Tiling CRISPR" to identify lncRNA functional domains. Using this approach, we identified Xist A-Repeats as the silencing domain, an observation in agreement with published work, suggesting Tiling CRISPR feasibility. Mechanistic analysis suggested a novel function for Xist A-repeats in promoting Xist transcription. Overall, our method allows mapping of lncRNA functional domains in an unbiased and potentially high-throughput manner to facilitate the understanding of lncRNA functions.

Drosophila Sister-of-Sex-lethal reinforces a male-specific gene expression pattern by controlling Sex-lethal alternative splicing.

In Drosophila, female development is governed by a single RNA-binding protein, Sex-lethal (Sxl), that controls the expression of key factors involved in dosage compensation, germline homeostasis and the establishment of female morphology and behaviour. Sxl expression in female flies is maintained by an auto-regulatory, positive feedback loop with Sxl controlling splicing of its own mRNA. Until now, it remained unclear how males prevent accidental triggering of the Sxl expression cascade and protect themselves against runaway protein production. Here, we identify the protein Sister-of-Sex-lethal (Ssx) as an inhibitor of Sxl auto-regulatory splicing. Sxl and Ssx have a comparable RNA-binding specificity and compete for binding to RNA regulatory elements present in the Sxl transcript. In cultured Drosophila cells, Sxl-induced changes to alternative splicing can be reverted by the expression of Ssx. Moreover, in adult male flies ablation of the ssx gene results in a low level of productive Sxl mRNA splicing and Sxl protein production in isolated, clonal cell populations. In sum, this demonstrates that Ssx safeguards male animals against Sxl protein production to reinforce a stable, male-specific gene expression pattern.

RNA element discovery from germ cell to blastocyst.

Recent studies have shown that tissue-specific transcriptomes contain multiple types of RNAs that are transcribed from intronic and intergenic sequences. The current study presents a tool for the discovery of transcribed, unannotated sequence elements from RNA-seq libraries. This RNA Element (RE) discovery algorithm (REDa) was applied to a spectrum of tissues and cells representing germline, embryonic, and somatic tissues and examined as a function of differentiation through the first set of cell divisions of human development. This highlighted extensive transcription throughout the genome, yielding previously unidentified human spermatogenic RNAs. Both exonic and novel X-chromosome REs were subject to robust meiotic sex chromosome inactivation, although an extensive de-repression occurred in the post-meiotic stages of spermatogenesis. Surprisingly, 2.4% of the 10,395 X chromosome exonic REs were present in mature sperm. Transcribed genomic repetitive sequences, including simple centromeric repeats, HERVE and HSAT1, were also shown to be associated with RE expression during spermatogenesis. These results suggest that pervasive intergenic repetitive sequence expression during human spermatogenesis may play a role in regulating chromatin dynamics. Repetitive REs switching repeat classes during differentiation upon fertilization and embryonic genome activation was evident.

Insertion position as well as the inserted TRS and gene sequences differentially affect the retention of foreign gene expression by simian hemorrhagic fever virus (SHFV).

Recombinant SHFV infectious cDNA clones expressing a foreign gene from an additional sg mRNA were constructed. Two 3' genomic region sites, between ORF4' and ORF2b and between ORF4 and ORF5, were utilized for insertion of the myxoma M013 gene with a C-terminal V5 tag followed by one of the three inserted transcription regulatory sequences (TRS), TRS2', TRS4' or TRS7. M013 insertion at the ORF4'/ORF2b site but not the ORF4/ORF5 site generated progeny virus but only the recombinant virus with an inserted TRS2' retained the entire M013 gene through passage four. Insertion of an auto-fluorescent protein gene, iLOV, with an inserted TRS2' at the ORF4'/ORF2b site, generated viable progeny virus. iLOV expression was maintained through passage eight. Although regulation of SHFV subgenomic RNA synthesis is complex, the ORF4'/ORF2b site, which is located between the two sets of minor structural proteins, is able to tolerate foreign gene insertion.

The conserved regulatory RNA RsaE down-regulates the arginine degradation pathway in Staphylococcus aureus.

RsaE is a regulatory RNA highly conserved amongst Firmicutes that lowers the amount of mRNAs associated with the TCA cycle and folate metabolism. A search for new RsaE targets in Staphylococcus aureus revealed that in addition to previously described substrates, RsaE down-regulates several genes associated with arginine catabolism. In particular, RsaE targets the arginase rocF mRNA via direct interactions involving G-rich motifs. Two duplicated C-rich motifs of RsaE can independently downregulate rocF expression. The faster growth rate of ΔrsaE compared to its parental strain in media containing amino acids as sole carbon source points to an underlying role for RsaE in amino acid catabolism. Collectively, the data support a model in which RsaE acts as a global regulator of functions associated with metabolic adaptation.